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Abstract 

In this work, we calculate the current distribution, in the close vicinity of the quan- 
tum point contacts (QPCs), taking into account the Coulomb interaction. In the 
first step, we calculate the bare confinement potential of a generic QPC and, in the 
presence of a perpendicular magnetic field, obtain the positions of the incompress- 
ible edge states (lES) taking into account electron-electron interaction within the 
Thomas-Fermi theory of screening. Using a local version of the Ohm's law, together 
with a relevant conductivity model, we also calculate the current distribution. We 
observe that, the imposed external current is confined locally into the incompress- 
ible strips. Our calculations demonstrate that, the inclusion of the electron-electron 
interaction, strongly changes the general picture of the transport through the QPCs. 
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1 Introduction 



At low temperatures, low-dimensional electron systems manifest peculiar quantum- 
transport properties. One of the key elements of such transport systems are 
the quantum point contacts (QPC) constructed on a two-dimensional elec- 
tron system (2DES). The wide variety of the experiments concerning QPCs 
[1], including quantum Hall effect (QHE) based Mach-Zehnder interferome- 
ter (MZI) [2,3], have attracted many theoreticians to investigate their elec- 
trostatic [4] and transport properties [5,6]. However, a realistic modelling of 
QPCs that also takes into account the involved interaction effects is still un- 
der debate. The magneto-transport properties of such narrow constrictions 
is typically based on the standard IDES [7], which relates the conductance 
through the structure to its scattering characteristics, considering typically 
a hard-wall confinement potential. The reliability of such non-interacting ap- 
proaches is limited, since interactions are inevitable in many cases and plays a 
major role in determining the electronic and transport properties. In order to 
account for the interactions simplified models are used with some phenomeno- 
logical parameters, which is not always evident whether such a description is 
sufficient to reproduce the essential physics. 

The QHE based MZI [2] has become a central interest to the community, since 
it provides the possibility to infer interaction mechanisms and dephasing [8,9] 
between the ES by achieving extreme contrast interference oscillations. In 
these experiments ESs [7,10,11,12] are assumed to behave like optical beams, 
whereas QPCs simulate the semi-transparent mirror in its optical counter- 
part. The unexpected behavior of interfering electrons, such as path- length- 
independent visibility oscillations, is believed to be related to long range e — e 
interactions. Thus, the experimental findings present a clear demonstration 
of the breakdown of, commonly used, Landauer's conductance picture away 
from the linear regime. Here, we calculate the effective potential in a self- 
consistent manner and , in addition, using a local version of the Ohm's law 
within and out-of-the-linear-response regime, we obtain the current distribu- 
tion near the QPCs. We essentially show that, in the presence of an lES, the 
imposed current is confined to this region, otherwise is distributed classically. 



2 Model, Results and Discussion 

Our aim is to calculate the distribution of the ES within a interacting model. 
We start with the bare confinement potential obtained from the lithographi- 
cally defined construction, following Ref's [4,13]. For a given pattern of (metal- 
lic) gates residing on the surface and the potential values Vg{x,y,0), one can 
obtain the potential experienced Vext{x,y) by the 2DES beneath using semi- 



2 



analytical scheme [13] yielding 



where k is the dielectric constant of the hetero-structure 12.4 for GaAs/AlGaAs) 
and r = {x,y). Given the external potential in the position space, it is 
straight forward to calculate the screened potential in the momentum space 
{q) by Vscr{q) = Vcxt{q)/^{q) using the Thomas-Fermi dielectric function, 
e{q) = 1 + l/{ao\q\), where ao/2 = = Kh^ /{me^) (for GaAs = 9.8 nm). 
We use this potential to initiahze the self-consistent scheme described below, 
to obtain density nci(r) and potential distribution in the presence of a perpen- 
dicular magnetic field, 5, at a finite temperature T, /U*(r) is SC'ly found. In 
the absence of an fixed external current /, /U*(r) is position independent and is 
constant all over the sample, which is in turn determined by the average elec- 
tron (surface) number density, riei(r). In our calculations we set riei = 3.0- 10^^ 
cm~^, corresponding to a Fermi energy E-p ~ 10.7 meV. Starting from V^cr(r) 
one can obtain the electron density distribution, within the Thomas-Fermi 
approximation [14] (TFA), from 



(1) 



nei(r) = j dED{E)f{[E + V{r) - fx^yksT), (2) 

where D{E) is the Gaussian broadened (single-particle) density of states (DOS) 
and /(^) = [l-|-e^]~^ the Fermi function. The total potential energy of an elec- 
tron, V{r) — Kxt(r) + Vff(i"), differs from the Hartree potential energy Vff(r) 
by the contribution due to external potentials and is calculated from 



y(r) = V;xt(r) + ^ / dr'K{r, r') nei(rO. (3) 

K J A 

For periodic boundary conditions, the kernel K{r, r') can be found in a well 
known text book [15], otherwise has to be solved numerically. 

In a classical manner, if a current is driven in y direction in the presence of a 
perpendicular B field a Hall potential develops in the x direction. Therefore 
the electrochemical potential has to be modified due to external field in the 
y direction with E{r) = V/i*(r)/e = /3(r)j(r), for a given resistivity tensor 
and boundary conditions, in the thermal-equilibrium (locally), which brings 
a new self-consistent loop to our problem. We calculate the electric field by 
solving the equation of continuity under static conditions, V.j(r) = 0, and 
V X E{r) — 0, for a fixed total current, self-consistently. 

As mentioned above, in the single particle model it is believed that the current 
is carried by the baUistic ID Landauer-Biittiker (LB)-ES and the conductance 
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is obtained by the transmission coefficients. Here we calculate, the electron and 
current density self-consistently and observe the different distributions of the 
lES under quantum Hall conditions, within the linear response regime. We 
will always consider the case, where Ep is larger than the hight of the barrier 
at the center of the QPC, so that the conductance is finite and it is at least 
equal to the first Landau energy {Q/2 = huJc/2) = eB/2m). In Fig. 1, we 
selected three representative B field values, such that (i) the system is almost 
compressible (a); (ii) the lES merge at the opening of the QPC (b) and; (iii) 
the lES percolate through the constraint (c). From the "classical" current 
point of view we observe that, the current biased from bottom is (almost) 
homogeneously distributed all over the sample if there exists no lES inside 
the constraint, the current passes through the QPC and ends at the right top 
contact, mostly (d). Fig. le, shows us that, the current is confined to the lES 
and conductance is quantized, whereas for Q = 1.4 it is a bit larger than e'^/h. 
These results indicate that, the ESs present structures inside the QPCs if one 
models them in a more realistic scheme rather than as a single point, although 
the conductance quantization remains unaffected. Since now we can calculate 
the widths of the lES, depending on the magnetic field and sample structure 
it is also possible within this model to obtain the electron velocity inside lES, 
which may be combined with a recent work by 1. Neder [16]. This work is 
based on non-Gaussian noise measurements at the Mach-Zehnder interference 
experiments. Their main finding is that the unexpected visibility oscillations 
observed, can be explained by the interaction between the "detector" and 
"interference" edge channels. The essential parameters of this work are the 
electron velocity at the detector edge channel and the coupling (interaction) 
strength between the detector and the interference channels. We believe that, 
an extension of our present model to a realistic system may contribute to the 
understanding of the mentioned experiments. 

We would like to thank R. R. Gerhardts for his fruitful lectures on screening 
theory, enabling us to understand the basics. The authors acknowledge the 
support of the Marmaris Institute of Theoretical and Applied Physics (ITAP), 
TUBITAK grant 105T110, TUBAP-739-754-759, SFB631 and DIP. 
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Fig. 1. The color coded iy{x,y) for = 1.4,1.24,1.13 (a-c) and correspond- 
ing current distribution indicated by the arrows. The calculations are done 
at Q,/kBT !v 0.025 for a fixed external current driven in the y direction 
(j(r) = jo(0, 0.1)), where the ID current density is set to —0.42 x 10~^ A/m. The 
gates defining the QPC, arc taken to be 300 nm apart and biased with —0.3 V. 
The 2DES is 85 nm below the surface and unit cell is a square with dimensions 
1500 X 1500 nm2. 
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